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We propose a framework to construct the ground-state energy and density matrix of an A'^-electron 
system by solving selfconsistently a set of single-particle equations. The method can be viewed as a 
non-trivial extension of the Kohn-Sham scheme (which is embedded as a special case). It is based on 
separating the Green's function into a quasi-particle part and a background part, and expressing only 
the background part as a functional of the density matrix. The calculated single-particle energies 
^ I' and wave functions have a clear physical interpretation as quasiparticle energies and orbitals. 
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I. INTRODUCTION 



The power of the Kohn-Sham (KS) implementation JJ of density functional theory (DFT lies in its ability to 
incorporate complex many-body correlations (beyond Hartree-Fock) in a computational framework that is not any 
^ , more difficult than the Hartree-Fock (HF) equations. In practice this means one has to solve single-particle Schrodinger 
' equations, with local or non-local potentials, in an iterative self-consistency loop. This simplicity makes KS-DFT the 
O ' only feasible approach in many modern applications of electronic structure theory. There is therefore continuing 
. interest, not only in developing new and more accurate functionals, but also in studying conceptual improvements 
y ' and extensions to the DFT framework. 

C/2 The present implementations of DFT can handle short-range interelectronic correlations quite well, but often 

fail when dealing with near-degenerate systems characterized by a small particle-hole gap. This seems to indicate 
that KS-DFT does not accurately describe the Fermi surface if it deviates significantly from the noninteracting 
one Umillllllll. In this respect, one of the most glaring inadequacies of KS-DFT is the fact that the 
physically important concept of quasiparticles is missing. Even formal knowledge of the exact exchange-correlation 
energy functional would only lead to the total energy and the density, since the individual KS orbitals have no special 
significance. An exception is the HOMO orbital and energy which govern the asymptotic tail of the density. It has also 
been observed numerically that in exact KS-DFT the occupied single-particle energies resemble the exact ionization 
energies also for the deeper valence hole states [lOl ; the deviations have been analyzed theoretically in terms of a 
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! reponse contribution to the KS potential 0, Il2j|. 
' Quasiparticle (QP) excitations in the Landau-Migdal sense form a well-known concept in many-body physics. 
' They are most readily understood as the relics of the single-particle (s.p.) excitations in the noninteracting system 
■ when the interaction is turned on (l5ifl6i|. In most electronic systems, or more generally in all normal Fermi systems, 
c/3 ' the bulk of the s.p. strength (i.e. the transition strength related to the removal or addition of a particle) is concentrated 
• in QP states. Especially near the Fermi surface, the QP states represent the dominant physical feature and should 
c/3 be described properly in any appropriate single-particle theory. 

The QP orbitals form a complete, linear independent, but generally nonorthogonal set. The completeness and 
linear independence follows from the fact that the QP states evolve from the set of single-particle eigenstates of a 
noninteracting Hamiltonian. Near the Fermi surface they coincide with the electron attachment states or with the 
dominant ionization states. Further away from the Fermi surface, QP states may acquire a width and correspond 
not to a single state but to a group of states in the (TV ± l)-electron system characterized by rather pure one-hole or 
one-particle structure. The QP states have reduced s.p. strength (i.e. the normalization of the QP orbitals is less than 



■ unity). Note that QP orbitals and str eng ths (at least near the Fermi surface) are experimentally accessible using, 
e.g., electron momentum spectroscopy \IT\ - 

It is the aim of the present paper to develop formally exact single-particle equations whose solutions can be 
interpreted as the QP energies and orbitals, and to explore how approximations can be introduced by the modelling 
of small quantities in terms of functionals. The resulting formalism will be called QP-DFT and yields, apart from the 
QP orbitals and energies, also the total energy and the density matrix of the system. 

To the best of our knowledge this is the first time such a formalism, aiming directly at the quasiparticle properties, 
is proposed. There are some similarities to the "correlated one-particle theory" as developed recently by Beste et 
ai. [n m, in that an energy-independent single-particle Hamiltonian is obtained whose eigenvalues are related 
to ionization energies and electron affinities. The interesting proposal in Ref. 0, is to construct such a s.p. 
Hamiltonian by means of a Fock-space coupled cluster expansion, separately in the N — 1 and N + 1 electron system, 
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and supplemented by a maximal overlap condition between the reference Slater determinant and the exact N electron 
ground state (i.e. Brueckner coupled cluster). However, it is not obvious how the overcomplete set of Dyson orbitals 
for the iV — 1 electron system is reduced to N single-particle orbitals because the technique seems to require that in 
general precisely N exact eigenstates of the iV — 1 electron system are projected out as principal ionization states. 
This may be problematic in cases where the s.p. strength is fragmented, or when considering the electron gas limit. 
In our opinion, a general formulation should be expressed using the QP states and their associated width. 

This paper is organized as follows. Sec. ITU reviews all necessary theoretical ingredients of the QP-DFT formalism. 
The natural framework for discussing QP properties is the propagator or Green's function (GF) formulation of 
mam (^odv theory [H A summary of some GF results is therefore provided in Sees. Ill AIII Bl See, e.g., 

[2^. '22I I23L l23 | for reviews and general papers on the specific use of GF theory in electronic structure problems. In 
Sec. Ill Gi the standard quasiparticle concept is introduced. The fundamental equations of the new QP-DFT formalism 
are established in Sec. IIII Al A discussion of the main properties of QP-DFT is then provided in the remainder of 
Sec. mil i.e. the embedding of both HF and KS-DFT, the potential exactness of the method, the asymptotic properties 
in coordinate space, and the electron gas limit. The final Sec. IIVI contains further discussion of QP-DFT (possible 
implementations, advantages and drawbacks), and concludes with a summary of the paper. 



II. GREEN'S FUNCTION THEORY AND THE QUASIPARTICLE CONCEPT 

In this section some well-established results and concepts are reviewed, in order to introduce notation and to provide 
an overview of the expressions needed for the QP-DFT formalism. The basic relations in GF theory can be found in 
many textbooks, e.g. [nlfl^- In Sec. Ill Gl the quasiparticle concept as introduced in is discussed; this material 

is also extensively treated in [l^ I25I |26|. 



A. Single-particle propagator 

We initially keep the discussion as general as possible, and consider a normal (non-superconducting) Fermi system 
with Hamiltonian H = Hq + V, where Hq contains the kinetic energy and external potential and V is the two-particle 
interaction. The s.p. propagator in the energy representation is defined as 



G{a,f3;E) = {^^\a^ . al + a}, . a^K) (1) 

E ~{H - E^) + 17] ' ' E+{H - E^) - ir] 



where a,f3, .. label the elements of a complete orthonormal basis set of s.p. states, the second-quantization operators 
tta (ajj) remove (add) a particle in state a (/3), and 7y > is an infinitesimal convergence parameter. The exact ground 
state of the A^-particle system is denoted by l^*^) and its energy by Eq . 

Inserting complete sets of eigenstates of H into Eq. leads to the standard Lehmann representation of the s.p. 
propagator, 

' ^ E^{E!^+'~E^) + ^7] E+{E,^-'-E^)~^7] 

where the l^"^^^) are the eigenstates, and E^^^ the eigenenergies, in the {N ± l)-particle system. The notation in 
the last line of Eq. 10), i.e. ej^*' — ±{E^'^^ — Eq ) for the poles of the propagator, and 

(4+')a-«ka|*r'); (z<-'). = (*rV„|0 (3) 

for the s.p. transition amplitudes, will be used throughout the paper. Note that the amplitudes are usually called 
Dyson orbitals in the electronic context. 

The poles e<j+' and ejj"' of the propagator are located in the addition domain (eg^', -l-oo) and the removal domain 
(— 00,60"'), respectively. In a finite system both domains are separated by an energy interval (eQ"',eg^'). The width 
of the interval is the particle-hole gap, 

4"'-4-'=<+'-2i?o^ + i?o^-'>0, (4) 
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where positivity is guaranteed by the assumed convexity of the versus N curve. Only the interval is physically 
relevant, but for definiteness one can take the Fermi energy as the center of the interval, 

^F^\{et^+er) = \{E^+'-Er'). (5) 

In an infinite system one has Eq^' — ef,^' — ep. 

The (1-body) density matrix [A^'^'] and removal energy matrix [M*"'] can be expressed in terms of the propagator 

as 

/HE 
— e"'^G(a,/3;i?), 

n 

EE «|atK,7?]K)=^6<-'(z<-')„(z<-')^= f ^e^^^'EGia^p-E). (6) 

n 

The expressions in terms of the e*j'' and z^"' follow by inserting the complete set of {N — l)-particle eigenstates of H 
between the ajj and Oq, operators. The equalities on the right of Eq. © can then be obtained by contour integration 
in the complex .E-plane, where the e"'^ factor selects the pole contributions in the removal domain. 

Any one-body observable of interest can be calculated with the density matrix. The removal energy matrix allows 
in addition to calculate the total energy through the Migdal-Galitskii sum rule 

i?o^ = iTrace([i/o][A^<-'] + [M<-']), (7) 

which can be obtained by exploiting the fact that Trace[M<"'] — {^q\Ho + 2V\^^). The removal part of the 
propagator is sufficient for these purposes. However, only the (inverse of) the full propagator has a meaningful 
perturbative expansion, which takes the form of the Dyson equation, 

[G{E)]-^ = [G,{E)]-^ - (8) 

where [Gq{E)\ is the noninteracting propagator corresponding to the Hamiltonian Hq and [5](£')] is the (one fermion 
line) irreducible selfenergy. In an ab-initio calculation, the physics input is controlled by taking a suitable approxi- 
mation for the selfenergy, but here the reasoning is in terms of the exact selfenergy. The latter plays the role of an 
energy-dependent s.p. potential. For a discrete pole of the propagator, e.g., the amplitude obeys 

([ffo] + [S(4*')])4*'=e<*'e', (9) 

In Fig. ^ the general diagrammatic structure of the selfenergy is shown. The two distinct types also correspond 
respectively to the sum of all energy-independent, and all energy-dependent, contributions to the selfenergy. 

B. Spectral function 

In general, the single-particle spectral function is related to the propagator as 

[Sm = -i-sign(6^-£;)([G(£;)]-[G(£;)]t) (10) 

= E(4"')(4+')^^(i? - 4+') + Y.^z\;')^zi-y)U{E - e'-y). (ii) 

n n 

The zero'th and first energy-weighted moments obey the sum rules 

/ + 00 
dES{a,(3-E) = (^'^|{at,a,}|^'^) 
-OO 
/' + 00 

Mc«,/3 = / dEES{a,[3-E)^ {■^^\{a%[ac.,H]}\^^), (12) 




FIG. 1: General structure of the exact self-energy S(_B) in terms of the exact 4-point vertex function T. The wavy line 
represents the antisymmetrized interaction, the double directed line is the exact propagator. The first term is the sum of all 
static selfenergy contributions Vhf, the second term involving F groups all energy-dependent contributions. 

where the braces denote an anticommutator, e.g. {ap,aa} = Sa,p- The equahties on the right of Eq. H12|l follow 

[similar to Eq. ^] by inserting complete sets of eigenstates of the Hamiltonian H in both the (N + 1) and {N — 1) 
system, 

- E^""'(^""')(4"')^+E^""'(4-')(4-')^ (13) 

n n 

Note that the integrations in Eq. (|12|l are over the entire energy axis. If the integration is restricted to the removal 
domain (— cx),ei?) one retrieves the one-body density matrix [A^'^'J and removal energy matrix [M'"'] as defined in 
Eq. (jni, i.e. 

/ep p+ca 
dES{a,(3;E)+ dES{a, (3; E) ^ N'^j + N'^^^, (14) 

-OO J €f 

and similarly for Ma^p = M^y^ + 
Writing the Hamiltonian as 

as a^^^a\asau, (15) 

with 

{af3\V\j5)as = {af3\V\lS) - {af3\V\5^) (16) 

the antisymmetrized interaction matrix element, the (anti)commutator on the right of Eq. (|12|l can be worked out 
explicitly as 

{at , [a^,H]} - {a\Ho\P) + Y,{ai\Vmasa\as. (17) 

7(5 

As a result, it is possible to express the sumrules in Eq. (|12|l in closed form as 

Na,,f3 = <5„,0, or [N] = [/], (18) 
M„,^ = {a\Ho\l3) + Y,{o^l\VmasN',-\ or [M] = [Ha] + [VhfI (19) 

where [/] is the identity matrix. These expressions will be used extensively in the following. The second term, [Vf/F], 
in Eq. H19|l is the sum of all static (energy independent) selfenergy contributions, and has the form of the HE mean 
field, but evaluated with the exact density matrix [A^'^'j. A diagrammatical representation is provided by the first 
term in Eig. ^ 



C. Quasiparticles 



In normal Fermi systems, the bulk of the spectral strength is concentrated in quasiparticle (QP) states which, in the 
Landau-Migdal picture, evolve adiabatically from the A^± 1 eigenstates of Hq, and can be regarded as the elementary 
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s.p. excitations in the interacting system. In its simplest form the QP contribution to the propagator can be written 
as a modified noninteracting propagator, 

plE- EQj - iWQj ^. E - eQj + iWQj 

where WQj > characterizes the width of the QP excitation at energy eQj, and zqj is the corresponding QP orbital. 
The first term in Eq. I|20() corresponds to excitations in the {N — l)-particle system, as indicated by the location of the 
poles in the upper half- plane (Imi? > 0). It contains the N lowest QP energies egj {j = 1, .., N), for which egj < ep- 
The higher QP states {j — N + 1, .., cxo) in the second term of Eq. (|20|) correspond to the {N + 1) system, as indicated 
by the location of the poles in the lower half-plane {ImE < 0), and have eqj > ep- 

According to Eq. ({Tn|l , the spectral function corresponding to the QP propagator in Eq. (^0)) reads 

oo 

[Sq{E)] = Y.('Q,){ZQ,)^C^^^{E - eg,), (21) 
i=i 

where the (Breit-Wigner) distribution Ca{x) = ^^^^ is normalized to unity and has width A. 

The set of QP orbitals is complete and linearly independent (in view of the 1-1 correspondence to the noninteracting 
system), but in general not orthogonal. The QP orbitals are different from the HF ones because of the inclusion of 
selfenergy contributions beyond HE (see second term in Eig.nj. This energy dependent part of the selfenergy is also 
responsible for a reduction of the spectroscopic strength, Zq^zqj < 1. 

For infinite and homogenuous systems (e.g. the electron gas) the situation is particularly transparant, since all 
matrix quantities are diagonal in a plane-wave basis, and reduce to continuous functions of s.p. momentum p. For 
momenta close to the Fermi momentum (p ~ pp) the QP energy is the solution of eqip) = p^ / (2m) -|- ReE(p; eQ(p)), 
and coincides with the Fermi energy at pp, i.e. eQ{pp) — ep. The spectral function S{p; E) is dominated by a sharp 
QP peak at eqip), having a width wq{p). The width wq{p) is proportional to ImE(p; eglp)) and vanishes for p pp. 
While the QP excitations are only unambiguously defined near the Fermi surface, the concept can be extended to all 
momenta. 

Similar considerations also apply to finite systems, where the QP excitation corresponds to a state or group of states 
having rather pure s.p. character. The {N ± l)-particle states near the Fermi energy (where no poles of the selfenergy 
are present) are usually characterized by spectroscopic factors of the order of (but smaller than) unity. In closed-shell 
atoms, e.g., the first ionization state V&o ~^ (corresponding to the first pole in the removal domain of the propagator) 
typically has a spectroscopic factor of about 95% |23|. In that case the QP orbital in Eq. 1)20(1 coincides with the 
Dyson orbital. More deeply bound orbitals may acquire a width, but have similar summed strength concentrated 
near an average QP energy. The QP orbital in Eq. (|20|l is then representative for the transition amplitudes to these 
states. 

The QP contribution to the sumrules in Eq. H12|) is given by 

oo oo 

[^q] = E^Qj4j; [^q] = E^Q^-^Q^-4r (22) 
J=l J=l 

Note that the QP width does not contribute to the 0^^ and 1^^ moment, and drops out from all subsequent consid- 
erations. 

Finally, it should be noted that the simple ansatz in Eq. (|20|l does not provide a valid description of a QP excitation 
j far away from the QP energy eqj. The energy-dependence in Eq. (|21|l can then be slightly misleading, as the tails 
of the Breit-Wigner distribution extend beyond the Fermi energy. In reality the strength of a QP excitation either 
belongs to the iV — 1 or to the iV-f 1 system. The QP contribution to the separate [A^<*>] and [M'='='] should therefore 
not be determined directly by Eqs. ()14I21|1 . but rather by the location of the poles in the lower or upper halfplane, 
i.e. 

N 



r HF ^ 
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and similarly for 

oo oo 

[^q']= E 'qAp [^^q']= E ^Q^^qAj- (24) 

i=N+l j=N+l 

III. QUASIPARTICLE EQUATIONS 

A. Derivation of the QP-DFT equations 

For the following it is important to realize that, given arbitrary hermitian matrices [Nq\ and [Mq\ with [Nq\ 
positive-definite, one can always write the unique decomposition of Eq. H22|l. This can be achieved by constructing 
the unique basis that solves the (generalized) eigenvalue problem 

[Mq\uj = \ANq\uj- u][NQ\uk = <5j,fc, (25) 

where [Nq\ plays the role of a metric matrix; the QP energies and orbitals given by eqj — Xj and ZQj ~ [NQ]uj then 
fulfill Eq. 

The eigenvalue problem in Eq. (^5)1 can be considered as a set of s.p. equations determining the QP orbitals 
and energies. We now rewrite the unknown operators [Nq] and [Mq] in a more useful form that suggests possible 
approximation schemes. 

Since the QP contribution to the spectral function is dominant, it makes sense to isolate it, 

[SiE)] = [SQiE)] + [SB{E)], (26) 

and concentrate on the residual small 'background' contribution [Sb{E)]. In fact, the full energy dependence of 
[Sb{E)] is not needed, since one can apply the quasiparticle-background separation of Eq. H26|l to the zero'th and 
first energy-weighted moments as well, i.e. one has 

[N] ^ [Nq] + [Nb]; [M] = [Mq] + [A/b], (27) 

where the total energy integrals [Nb] — [Ng^] + [Ng^] and [Mb] = [-^-^s *] + [-^^s^'] can again be split in a removal 
and addition part. Note that the matrices [N] and [AI] on the left side of Eq. H27|l are known in closed form through 
Eqs. (|18I19|) . so it follows that 

[Nq] = [I]-[Nb], (28) 
[Mq] = [Ho] + [Vhf]~[Mb]. (29) 

One then arrives at the remarkable conclusion that modelling the background contributions [M]j^'] and [A'^^j^'] as 
a functional of e.g. the density matrix [A^*"'], is sufficient to generate a selfconsistent set of s.p. equations. Using 
Eqs. (|28I29() the eigenvalue problem in Eq. (|25(l can be expressed as 

{[Ho] + [Vhf{N'-^}] [Mb{N'-^}])u, = A,([/] - [Nb{N^-^}])u,, (30) 

where the functional dependency of [M b] and [Nb] is indicated between braces. Note that also the HF-like potential 
[Vhf] is by definition [see Eq. (|19(l ] expressed in terms of the density matrix, as indicated in Eq. (|30|l . 

Having an initial estimate for [TV'"'] allows to construct the matrices [Nb] and [Ms], as well as the HF-like 
potential [Vhf]- The eigenvalue problem in Eq. H3UI) can now be solved, yielding QP energies eQj — Xj and orbitals 
-^Qj ~ (W ~ • The N solutions with lowest energy represent excitations in the iV — 1 system, and 

should be used to update the density matrix [TV*"'] = [Nq^] + [N^^], 

N 

[A^U] = E^Q^-4 + iN^B-'iN'-'}]- (31) 
This closes the selfconsistency loop, which can be iterated to convergence. The total energy then follows from Eq. 0, 



1 ^ 

< = 2 ^ ^-^I^"] + "^'^'"^^ + i^'^"" iiHo][K-'{N'~'}] + [A4-'{^'-'}]) • (32) 
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The above formalism, henceforth called quasiparticle DFT (QP-DFT), generates the total energy, the density 
matrix, and the individual QP energies and orbitals, starting from a model for the background contributions [-M^"^'] 
and [A''^*] as a functional of the density matrix. It is intuitively clear that this is a reasonable strategy: the external 
potential appears directly in the QP hamiltonian [M] — [Mb] through the [Hq] term of Eq. (|19|l . and primarily 
influences the position of the QP peaks. One may then assume the background part to be generated by 'universal' 
electron-electron correlations, and to be a good candidate for modelling. Equation determining the QP energies 
and orbitals, is the central result of the present paper, and its properties will now be discussed in detail. 



B. HF and KS-DFT as special cases of QP-DFT 



We first show that both HF and KS-DFT theory are included in the general QP-DFT treatment. 

For HF this is rather obviously achieved by setting all background quantities [M^^'] and [N'j^^] equal to zero in 
Eq. H30|l . Since the metric matrix on the right of Eq. (|30|l is now simply the identity matrix, one has zqj = Uj and 
the QP orbitals form an orthonormal set obeying {[Hq] + [VHF])zQj = egjZQj. The matrix [Vhf] is given by Eq. itT^ 
where the density matrix [A'^*"'] in the present approximation follows from Eq. (|31|l with [-/V^ '] = 0, i.e. 

N 

[^'-'] = E^Q^-4- (33) 

One can see that [AI] = [Hq] + [Vhf] is just the ordinary HF hamiltonian, and the total energy obtained from Eq. H32f) 
by setting [7V<j-'] = [M^g^] = 0, 



^lj2^^Qj[Ho]zQ,+eQ,) (34) 

is equivalent to the HF total energy. 

The KS-DFT case is somewhat more difficult, but one also starts by setting [iV^"^'] = in Eq. H3U|I . so the QP states 
ZQj form an orthonormal set obeying 

{[Ho] + [V„] + [Vf] ~ [Mb])zq, = CQ.ZQ,. (35) 

Here the HF-like potential has been split into its direct and exchange components, [Vhf] = [Vh] + [Vf]- As the 
density matrix is again given by Eq. (|33|l . both components only receive contributions from the occupied {j = 1, .., A'^) 
orbitals zqj: 

N 

{VhU^ = Y,{aj[vmY,i^QM^Qj);, 

7(5 j=l 

N 

{Vfu,^ = -Y.{aj[vmY,{^QM^Qj);- m 

7(5 j=l 

For the Coulomb interaction and taking as s.p. labels for the electrons the space coordinate and (third component 
of) spin, a = rrus , this reduces to the familiar expression 

Note that for compactness we continue to employ the general matrix notation used so far, with the understanding 
that sums over s.p. labels should be replaced by coordinate space integrations where appropriate. The total energy 
follows from Eq. ^ with [N^g^] = 0, 



1 ^ 1 

Eo = 2 E(4 [^ol^Qj- + ^Q^) + ;jTrace[M<j-'] . (38) 



The unknown [Af^j"^'] should now be determined by identification with the results of a KS-DFT calculation. 
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We allow explicit dependence on the KS orbitals (flxsj, and write the exchange-correlation energy functional as 

N 

Exc = "^KSj [<^xc]<fKSj- (39) 

i=i 

Specializing e.g. to a hybrid functional (the derivation for a functional with general orbital dependencies proceeds in 
a similar fashion) with a fraction f3 of exact exchange, the matrix [exc] in coordinate-spin space reads, 

<Pi<'5i(rm-s)<j5't-o,-(r'TO') 

e,,(r m„ r' m^) = -f ^ |r + '5,„^.„,, 5(r - r')/Mp(r), Vp(r)), (40) 

where the second term contains a local functional of the electron density p(r) and its gradient (assuming for simplicity 
a spin-saturated system). The corresponding exchange-correlation potential [Vxc], appearing in the KS equation, then 
reads 

Identification of the orbitals and energies in the KS equation 

{[Ho] + [Vh] + [VxcDfRSj = eKSjfKSj (42) 

with those in Eq. then requires 

[Mb] = [Vf] - [Vxc]. (43) 

Identification of the KS total energy 

N 



Eo = E ^KsA^o] + [Vh])^ks, + Exc 

1 ^ 1 ^ 

= ^Yi'PKSjiHo]'PKSj + <^KS]) + Exc- V^KSj['^xc]^KS] (44) 



with Eq. requires 



i=i j=i 



N 



Trace[M<j->] = ^ V'L5j(2[e.c] - [Vxc])^kSj. (45) 



Choosing the [Mg ] operators as 



[M<,-'] = [P]{2[exc] [Vxc])[P]; [M'g^] = [Vp] [Vxc] [P](2[e..] - [Vxc])[P], (46) 

where [P] = X^jLi ^Qj'^Qj projects onto the occupied QP orbitals, fulfills the requirements of both Eq. (|43|) and of 
Eq. 1)45(1 ; this choice therefore leads to the same results as the KS-DFT approximation, for the total energy as well as 
the orbitals and orbital energies. 

One concludes that the QP-DFT formulation is flexible enough to reproduce HE or KS-DFT results by speciflc 
choices of and [Afjj"^']. Note that the KS-DET formalism is embedded in a slightly tortuous way, as [A^g*'] = 

would normally imply the absence of a background contribution to the first-order moment [Mj^"^'] as well. A similar 
situation is encountered in Brueckner-Hartree-Fock theory j2^, in which the real part of the selfenergy causes a shift in 
the QP energy, but the corresponding reduction of strength through the imaginary part of the selfenergy is neglected. 



C. Potential exactness 



It should be clear that GF quantities like the propagator, spectral function, or selfenergy, are well-defined and can 
in principle be calculated exactly. Also the separation of the spectral strength into QP and background parts can be 
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performed for any system. In fact, defining the QP part can usually be done in several ways which are all, however, 
equivalent at the Fermi surface. Note that the density matrix and the total energy only contain the O^'^ and 1^* 
energy-weighted moment [A''*"'] and [M<~'], and do not depend on the QP-background separation. As long as one 
sticks to one unique prescription, the QP orbitals and energies are well-defined quantities, and so are the background 
contributions [N'^^] and [M^*']. The fundamental expressions in Eq. (|30I32|I are therefore always valid. 

As an example of an unambiguous definition of QP excitations one can e.g., in the spirit of Landau-Migdal theory, 
consider the lowest order expansion of the selfenergy around ep: 

[E(i?)] = [EM] + (i?-6f)[E'M], (47) 

where is the first derivative of the selfenergy with respect to energy. For infinite homogenuous systems the 

right side of Eq. H47|) is a real function, since lm'S{p, E) = ^Iml](p, = at the Fermi energy E = ep. In a finite 
system the right of Eq. 147(1 is an hermitian matrix, since = in a finite pole-free interval around the 

Fermi energy. Substituting the linearized selfenergy l(47|l into Eq. © leads to 

([i/o] + [S(eF)] - ^F[^'{ep)])u = E{[I] - [J:'{ep)])u, (48) 

defining the QP excitations. Note the formal similarity to Eq. H3U|) . and the fact that the metric matrix is positive- 
definite as [S'(ei?)] is a negative-definite hermitian matrix. While Eq. H48() results in a correct description for E « ep, 
it may not be appropriate away from the Fermi surface, when the linear approximation ((47|l breaks down. 

It should be stressed that the possibility of defining QP excitations in various ways is not a shortcoming of the 
present approach, but rather a reflection of the physical reality that QP excitations are only unambiguously defined 
near the Fermi surface, where the density of states in the N ±1 system is small. The complete description of the 
s.p. properties in an interacting many-body system is contained in the energy dependence of the spectral function. 
Any effective s.p. Hamiltonian can at most describe the peaks in the spectral function, i.e. identify the energy regions 
where the s.p. strength is concentrated and assign an average transition amplitude to this region. 

Finally we should explain the use of the density matrix as the independent variable for the functional modelling of 
[iV^"^'] and [M^*'] in Eq. 1(30(1 . In the present context it seems most natural to do this in terms of a matrix quantity. 
Moreover, for the class of systems with fixed two-body interaction V and varying nonlocal external potentials (i.e. 
varying Hq), the fundamental theorem of density matrix functional theory p9l | implies that all quantities can be 
expressed as functionals of the density matrix. This holds in particular for the [A'^^j^'] and [il/^"^'], once a unique 
prescription is agreed upon. 

It is of course also possible to consider the electron density p as the basic variable, provided one restricts oneself 
to the class of systems with local external potentials. In that case one is again formally assured of the existence of 
functionals for all quantities. Note that the functional in Sec. IIII 51 reproducing the KS-DFT results, even when the 
exact KS functional is used, is not the exact QP-DFT functional since it does not yield the correct density matrix; 
the exact QP-DFT fimctional would definitely have [N'j^''] ^ 0. 



D. Asymptotics in coordinate space for finite systems 



For notational convenience the spin dependence is dropped here. In coordinate space the background quanti- 
ties A''^~'(r,r') and M]j~-'(r, r') represent nonlocal operators, as they have a finite trace. The complete background 
operators A''B(r,r') and MB(r,r') can be local or nonlocal. 

It seems plausible that the background operators are short-ranged in the sense that correlation effects should become 
insignificant at large distances. I.e. the limit |r| — > oo of the background operators should drop out when considering 
the asymptotic behavior of the quasiparticles in Eq. 1(30(1 . This requires for |r| —^oo that 

I y"dr'iVs(r,r')z(r')|/|z(r)| ^0 (49) 

for an arbitrary s.p. orbital z(r). It is sufficient to take [Nb] local and asymptotically vanishing, but nonlocal 
implementations are also possible. For [Mb] one requires that for |r| — oo, 

AfB(r,r')-> ^^^ + i?B(r,r') (50) 

|r — r'l 

where Hb{v,y') is again short-ranged in the sense of Eq. ((49(1 . Note that one does not require [iV^~'] to be short- 
ranged, but its contribution to the Fock term in [VfrF{A^*"'}] of Eq. ((30(1 should be canceled by [M^] in the limit 
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|r| — > oo. The [Ng '] contribution to the Hartree term remains, since this is not affected when taking the hmit |r| oo 

in Eq. (|2DJ). 

With the characteristics H49I50|I of the exact QP-DFT functional, the background operators simply drop out of the 
asymptotic behavior of Eq. (|30|1 . and one is left with (for an external charge Z) 

/ 7 — AT \ ^ I" dr' 

\-2^-^--'Q^)'Q^^'^^T.^'^Q=^^) J ^^4.(r'K(r') = o. (51) 

Note that, for |r| — + oo, one has ZQi{r) — Ui{r) since the metric matrix [/] — [Nb] becomes the identity by virtue 
of Eq. H49|l . The Hartree term gives rise to the N/r potential, and only the Fock term built with the N lowest QP 
orbitals needs further analysis. In fact, one immediately sees that the behavior of the Fock term is the same as the 
one studied in Ref. (or more generally in Ref. 3(3|). The asymptotic expansion 

IT— = r + ^ ^TT^^cos..), (52) 

where Pl is a Legendre polynomial and uj the angle between r and r', can be substituted in the Fock term of Eq. H51|l. 
The first {L = 0) term in Eq. (|5^ then only contributes in case of removal-type (z = 1, ..,iV) QP orbitals, since [see 
Eq. ^] 

j dv' z^j{r')ui{r') = zj^jU, = ul[NQ]u^ = S^^, (53) 

and effectively changes the asymptotic charge to Z — + 1. A similar analysis as in Ref. [33| then proves that the 
removal QP orbital with the slowest decay is the HOMO one (z = N), and goes like ~ e~"''r'^ with = —2€qn and 
/3 = {Z — N + l)/a. The Fock term imposes the same exponential decay (except in very special cases) on the other 
removal QP orbitals, but with an accompanying power law that decreases faster by at least l/r^. 

For the addition-type {i = N + 1, .., +oo) QP orbitals the first term in Eq. 1)52(1 does not contribute to the Fock 
term in Eq. H51|l . and the asymptotic charge remains Z — N. This means that for neutral atoms no Rydberg sequence 
of bound addition-type QP orbitals is present, which is in agreement with the real physical situation. Also note that, 
unlike HE, the presence of the [Mb] term in Eq. (|30|l does allow for the possibility that some addition-type QP orbitals 
become bound for neutral systems. The asymptotics of the bound addition-type QP orbitals is again correct, with an 
exponential decay governed by their individual separation energies cq,. They go like ~ Q-ar^fS i^if^i^ q,2 _ —2eQi and 
P = {Z — N)/a; the Fock term in Eq. 1)51(1 cannot influence this, since all removal-type QP orbitals decay faster. 



E. The electron gas limit 



We consider for simplicity the spin-unpolarized electron gas at density p, and suppress the spin indices. To find 
the electron gas limit of the [M^"*"'] and [N^j^''] quantities one needs the momentum distribution n{p) and the removal 
energy r(p), defined as 

n(p) = (vE-ola^aplvE-o); r{p) ^ {^oWpK, H]\^o) ■ (54) 

The QP-background separation can be made once a QP spectrum eqij)) and strength Zq{p) is defined. Near the 
Fermi surface the standard definitions 

eQ(p) =pV(2m) +ReS(p;eQ(p)); Zq(p) = 1/(1 - ^^""f^'^^ E^e^ip) (55) 

are valid. For momenta far from pp Eq. ((55(1 may break down, e.g. the equation defining eqip) may have multiple 
roots, or the strength Zq{p) may become larger than unity. An alternative definition, valid for all momenta, must 
then be used. For an example of such a procedure we refer to |32| . where the QP-background separation was made in 
the context of a GW description [33Ll34| of the electron gas. 

The background quantities in the electron gas are now expressed as 

N'b\p) = n{p) - e{pF - p)Zq[p)- N'^\p) = 1 - n{p) - 0{p - Pf)Zq{p) 

M'g \p) = rip) - 9{pF - p)ZQ{p)eQ{p)- M'/>{p) = p^/{2m) + Vf{p) - r{p) - e{p - pF)ZQ{p)eQ{p), (56) 
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where Vpip) is the Fock-hke potential (evaluated with the exact momentum distribution n(p)). The typical behavior 
of these background quantities as a function of momentum can be found in Fig. 1 of Ref . 32] . 

For the purpose of applying the electron gas results to inhomogenuous systems, it is not a priori clear what the 
equivalent of the Local Density Approximation would be. Working with the Wigner transform would seem the most 
straightforward, 



Jp-(r-r') Ar(-) 



(2n) 



(57) 



p=p(R) 

but there are several possibilities which reduce to the same electron gas result in the homogenuous limit. 



IV. SUMMARY AND DISCUSSION 



The concept of quasiparticles is an important tool to understand and describe normal Fermi systems. In this paper 
we developed a set of single-particle equations whose solutions correspond to the QP orbitals and energies. When 
the residual small background contributions are expressed as universal functionals of the density or density matrix, a 
single-particle selfconsistency problem (the QP-DFT scheme) is generated that can be easily solved for an approximate 
choice of the functionals. 

The QP-DFT scheme would seem to offer many advantages as compared to KS-DFT. There is no need for the 
difference between the kinetic energy of the interacting systems and a reference system. All s.p. orbitals and energies 
have physical meaning, in contrast to the KS orbitals. The asymptotic behavior of the QP orbitals comes out correct, 
provided the background operators are short-ranged. On the down side: since we no longer have a sharp Fermi 
surface, particle-number conservation is not automatically guaranteed, and should be built into the functionals. 

The fact that KS-DFT is built in as a special case is a very important feature of QP-DFT. In a sense, one cannot 
do worse than KS-DFT, since one adds more parameters to the model. Moreover, the new parameters are truly new 
degrees of freedom (the introduction of the metric matrix, allowing a softening of the Fermi surface), which cannot 
be mimicked by taking a more sophisticated KS functional. 

The modelling of the background operators [N^ ^], [M^^] is basically virgin territory. One option would be to 
exploit the relation H46|l with the KS formalism, using an existing XC functional form for [M^'], adding a similar 
form for [A^^j"^'], and (re)parametrize by fitting to total energies and ionization energies in a training set of atoms and 
molecules. 

Alternatively, one could devise parametrizations by performing GF calculations on a series of test systems, and 
construct the background operators directly from the calculated spectral function. A step in this direction was taken 
in Ref. [s^, where we applied an ab-initio selfenergy of the GqWq type to a series of closed-shell atoms. The QP- 
DFT scheme outlined in the present paper was used in first iteration (no selfconsistency) to generate the first-order 
corrections to the HF picture. We then constructed a simple QP-DFT functional that depends only on HF quantities, 
but was able to reproduce the most important results of the underlying ab-initio model. 
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